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Abstract. - The antiferromagnetic molecular finite chain Cr6 was studied by inelastic neutron 
scattering. The observed magnetic excitations at 2.6 and 4.3 meV correspond, due to the open 
boundaries of a finite chain, to standing spin waves. The determined energy spectrum revealed an 
essentially classical spin structure. Hence, various spin-wave theories were investigated in order 
to assess their potential for describing the elementary excitations of finite spin systems. 
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Since the discovery of quantum tunneling of the mag- 
netization in the molecules Mni2 and Fes, the interest 
in molecular nanomagnets has soared [1]. Most of these 
molecules are characterized by dominant nearest-neighbor 
Heisenberg exchange interactions between the magnetic 
metal centers within a molecule, and additional weaker 
magnetic anisotropic terms (such as on-site anisotropy, 
dipole-dipole interactions, etc.). In these clusters, the na- 
ture of the ground state and elementary spin excitations 
is hence determined by the Heisenberg interactions. For 
a proper description of their magnetism it is thus of fun- 
damental importance to understand the general question, 
what spin ground states and excitations may emerge from 
Heisenberg interactions on finite spin lattices (where finite 
is to be understood as to mean a few ten spin centers at 
most). 

In this context, the antiferromagnetic rings and the Ke- 
plerate molecule Fe3o have become prototypical examples 
[2-6] . Their low- lying energy spectrum exhibits a remark- 
able structure: In an energy vs total spin S representation 
it consists of a set of rotational bands whose energies in- 
crease as E(S) cx S(S + 1) [7-10]. Physically, the lowest 
rotational band, the L band, reflects the rotational degree 
of freedom of the Neel-like antiferromagnetic ground state; 
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the higher rotational bands, collectively called E band, 
correspond to the (discrete) antiferromagnetic spin-wave 
excitations. For the case of the antiferromagnetic rings, 
this picture of the excitations has been confirmed exper- 
imentally in much detail [11], and also for F330 evidence 
is strong [12]. Interestingly, the rotational band structure 
is found in other finite spin lattices, too, and is generic 
in this sense. Its emergence is correlated to a basically 
classical spin structure [8, 10, 11]. 

So far, most of the understanding of the elementary 
excitations in Heisenberg spin clusters was inferred from 
exact numerical diagonalization studies, which are obvi- 
ously limited to systems with few spin sites. The classical 
spin structure identified in the above systems, however, 
suggests to consider approaches such as spin-wave theo- 
ries (SWTs), which are standard in 3D systems, and have 
been demonstrated to be very capable also in 2D and ID 
systems [13]. As a great advantage SWTs would allow 
one to treat much larger spin clusters, and exploring their 
applicability to molecular nanomagnets is hence of high 
interest, but no systematic study has emerged yet. Linear 
and modified SWTs were recently applied to Mni 2 and 
Fe3o [14-16], with no general conclusions. 

The examples of the antiferromagnetic rings, Fe3o, 
or M1112, exhibit relatively high symmetry and periodic 
boundary conditions. In this work, finite antiferromag- 
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Fig. 1: Magnetic susceptibility vs temperature, x(T), of Crg. 
Open dots: data; thick line: calculated x(T) as discussed in 
the text; thin lines: contributions from (a) the L-band states, 
(b) the remaining states, and (c) an impurity. The inset shows 
the molecular structure of the [Cr 6 Fu (0 2 CCMe3)io(H 2 0)] 3_ 
anion in Cr6. 

netic Heisenberg chains, as described by the spin Hamil- 
tonian 

N-l 

H = -J^Si-S i+1 , (1) 

are considered (J < is the coupling constant, TV is 
the number of spin centers with spin length Si). Fi- 
nite antiferromagnetic chains were studied extensively be- 
fore, for examples see Ref. [17]. Experimentally they 
were realized by intercepting infinite chains with non- 
magnetic dopants. This, however, yields a distribution 
of chain lengths, preventing observation of, e.g., individ- 
ual spin- wave excitations [17]. The molecular compound 
[ n Pr 2 NH2]3[Cr 6 F 11 (O 2 CCMe 3 ) 10 (H 2 O)] (or Cr 6 in short) 
[18], in contrast, epitomizes a perfect finite chain: The 
anion of Crg consists of six spin-3/2 Cr 3+ ions forming a 
horseshoe structure, see fig. 1 [19], and the generic finite- 
chain Hamiltonian eq. (1) is the obvious magnetic model 
(with N = 6, Si = 3/2). All molecules in a sample are 
identical, thus resolving the problem of the chain-length 
distribution, such that the Crg compound provides un- 
precedented opportunities for the experimental study of 
the excitations in finite antiferromagnetic chains. 

The spin excitations in finite chains show similarities 
to rings, but also differences. In particular, due to the 
open boundaries, the spin-wave excitations correspond to 
standing spin waves, as opposed to running waves in rings. 
In this work, we used inelastic neutron scattering (INS) to 
study the magnetic excitations in Crg experimentally, and 
hence to observe the standing spin waves in a finite antifer- 
romagnetic ID system for the first time. We also extended 
the rotational-band structure concept to finite antiferro- 
magnetic chains, and performed various SWT calculations 
to compare with the exact quantum mechanical result in 
order to explore the usefulness of these approaches. 
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Fig. 2: Neutron energy-loss data of Cr6 measured with (a) 
A = 3.8 A and (b) A = 3.0 A. In (a) the curves are shifted 
for clarity. The solid line represents the calculated spectrum 
at low temperatures (J = — 1.27meV, T = 1.6 K, linewidth 
= 0.2 meV). A± denotes the _E-band splitting, schematically 
depicted in the inset. 

A powder sample of undeuterated Cr 6 was synthesized 
as in [18]. The magnetic susceptibility was measured 
at a field of 0.1 T with a SQUID magnetometer (Quan- 
tum Design). INS experiments in the temperature range 
1.5 — 20 K were performed on the instruments IN5 (In- 
stitut Laue-Langevin, Grenoble, France), and IRIS (ISIS, 
Didcot, UK). The data treatment employed a correction 
for detector efficiency by a vanadium standard. The data 
shown in this work were obtained on IN5 with initial wave- 
lengths of A = 3.0 A and 3.8 A. The experimental resolu- 
tion at the elastic line was 0.40 meV and 0.17 meV, respec- 
tively. 

The magnetic susceptibility vs temperature, x(T), of 
Crg is shown in fig. 1. It exhibits a peak at 5K and 
a shoulder at 35 K. Observation of such a behavior is 
unusual; a similar behavior has been reported previously 
only for the antiferromagnetic wheel CrsNi [20]. The cal- 
culated susceptibility (thick solid line in fig. 1) was based 
on eq. (1) with J = — 1.27meV as found from INS (vide 
infra), g = 1.96 which is typical for Cr 3+ , and a Curie- 
type contribution which accounts for a small amount of 
paramagnetic impurity (C = 0.08 emu K mol -1 ). The 
agreement is very good, given that g and C were the only 
adjustable parameters. 

The INS data recorded on IN5 at various temperatures 
are shown in fig. 2. Six inelastic features at about 0.4 meV 
(I), 1.0 meV (ii), 1.7 meV (Ph), 1.9 meV (iii), 2.6 meV 
(IV), and 4.3 meV (V) are observed. The broad band Ph 
and the rising background are assigned to a non-magnetic 
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Fig. 3: Calculated energies vs total spin S for (a) Cr6 with 
J = — 1.27meV and (b) a N=6 spin-3/2 ring with J r i ng = § J- 
The lowest states for each S > form the L band (dark gray 
levels), the 4 next-higher levels for each S > 1 the E band 
(black levels), and the remaining states the quasi-continuum 
(light gray levels). 

origin (incoherent scattering of protons, phonons, inter- 
nal rotations of methyl groups, etc.) [21]. The other, 
resolution-limited peaks are assigned to magnetic transi- 
tions. The temperature dependence of the 3.8 A spectrum 
reveals that peaks I, IV, and V are cold transitions, while 
peaks ii and iii are hot transitions. 

The INS spectra were simulated by numerically diag- 
onalizing cq. (1) and calculating the intensities from the 
wavefunctions using the formulas of Ref. [22]. The solid 
line in fig. 2(a) represents the calculated spectrum for 
J = — 1.27 meV. The agreement with experiment is very 
good, considering that the model includes only one free 
parameter (J) [23]. The calculated position of peak V at 
4.1 meV is 5% too low, which indicates a weak variation 
of the exchange coupling constants along the chain [24]. 
These results demonstrate that the generic antiferromag- 
netic chain model eq. (1), while it does not describe Crg 
perfectly, clearly provides a very good starting point for 
the discussion of the physical nature of the elementary 
excitations in Crg. 

Figure 3(a) presents the calculated lowest energy levels 
vs S, with the observed transitions indicated by arrows. 
The three cold peaks I, IV and V correspond to transitions 
from the S = ground state to the five lowest triplets 
(which will be identified below as the S = 1 states of the 
L and E bands, respectively) . Peaks IV and V each consist 
of two transitions, which are too close to be resolved in the 
experiment. Peaks ii and iii originate from hot transitions 
within the L band, as indicated in fig. 3(a). 

In the following, the low-temperature triplet excitations 
I, IV, and V in Crg will be analyzed from two points of 



view. First, Cr@ is regarded as a hexanuclear antiferromag- 
netic ring with a missing bond ( Jei) as a perturbation, i.e, 
the Hamiltonian is written as TL = —J(Y]j_y Si ■ Si+i + 
Sq ■ Si) + V with V = J Si ■ Se- Second, various variants 
of SWT will be applied. 

The situation in an antiferromagnetic iV=6 spin-3/2 
ring is briefly recalled [10]. The energy spectrum for 
Jring = 5/6 J is depicted in fig. 3(b) (the factor 5/6 ad- 
justs for the different number of antiferromagnetic bonds 
in the chain and ring). It is characterized by a L band, 
an E band, and the quasi-continuum. The separation of 
the spectrum into L/E bands and a quasi-continuum orig- 
inates from the fact, that the INS intensity of transitions 
from states of the L band to states of the quasi-continuum 
is essentially zero, i.e., all intensity goes into transitions 
within the space of the L and E bands. Importantly, 
these properties are intimately connected to a classical 
spin structure [8, 10]. 

Comparing figs. 3(a) and 3(b) reveals that the concept 
of a L band, E band, and quasi-continuum is still use- 
ful for the finite chain. This is obvious for the L band, 
though not for the E band. However, a careful inspec- 
tion of the transition matrix elements shows that almost 
all low-temperature intensity goes into transitions to the 
states displayed in dark gray and black in fig. 3(a). For 
instance, transitions I, IV, and V account for 97.4% of the 
total intensity. Thus, although the states which form the 
E band in Crg are not easily recognized from the energy 
spectrum, they are clearly identified by the transition in- 
tensities. 

For a finite antiferromagnetic chain, as compared to an 
antiferromagnetic ring, two differences are noted. First, 
the gap between the ground state and the lowest triplet 
Ao (singlet-triplet gap) is rather small. This explains the 
peak and shoulder in x(T). Figure 1 disentangles the con- 
tributions to x{T) from (a) the L band (which is responsi- 
ble for the peak at about 5K), (b) all other states (which 
produces the shoulder at about 35 K), and (c) the impu- 
rity. As the L band is governed by Ao, the observation of 
a peak separated from a shoulder directly reflects a com- 
paratively smaller Ao: In the antiferromagnetic rings the 
L-band contribution (a) is shifted towards higher tempera- 
tures and hence less separated from contribution (b) , such 
that only a broad maximum is observed in x(^) P]- A 
second important difference is the splitting of the E band 
into two subgroups in the N=6 antiferromagnetic chain 
[from fig. 3(a) this is obvious for the triplets, but less so 
for S > 1]. The splitting of the E band is a direct con- 
sequence of the different boundary conditions for the spin 
waves in chains and rings. 

Considering a finite antiferromagnetic chain as a ring 
with a missing bond provides an intuitive picture. The 
cyclic symmetry of a ring implies a shift quantum number 
q = (r- f )^, with r = 1,2,..., N and f\q) = e~ tq \q) 
(T is the shift operator). The ground state belongs to 
q = 0, the triplet of the L band to q = 7T, and the triplets 
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Fig. 4: Spin-wave dispersions vs g obtained from linear SWT 
for chains with periodic (dotted line) and open (dashed line) 
boundary conditions. Symbols indicate the predicted energies 
for Cr6 (A) and a JV=6 ring (o). 



of the E band to q = ±571", ±§tt for a N = 6 ring. In 
the ring, the states with q and — q are degenerate by sym- 
metry, but in Ci'6 the perturbation of the missing bond 
produces a splitting in first order since (— q\V\ + q) ^ 0. 
Accordingly, the wavefunctions are no longer eigenstates 
of q (running waves), they are given by the symmetric and 
antisymmetric mixtures 



|g>± = | + g)±|- 9 ), 



(2) 



as depicted in the inset of fig. 2. These wavefunctions 
can be considered as the cosine and sine components built 
from I — q) and | + q), in analogy to e~ %q ± e lq , and hence 
represent standing waves. In Cr@ the splitting between 
them is given by A± = 2|(-f |V|§)| = 2|(-^|V|^)|, and 
the iS-band triplets organize into two groups separated by 
A±, see fig. 3a (the small splitting of the two states within 
one group, which increases with S, is a higher-order effect). 
Obviously, the observation of peaks IV and V in Crg, fig. 2, 
directly reflects the splitting A± due to the formation of 
the \q)± states, i.e., the observation of standing spin waves 
in a finite antiferromagnetic chain. 

Another simple rational for the -E-band splitting is 
obtained from linear SWT. The dispersion relations for 
an infinite antiferromagnetic chain with periodic or open 
boundary conditions (BC), respectively, are [7,25] 



periodic BC: 
open BC: 



3 ring sin(g)|, 
huo q = 2Si\Jsm(q/2)\. (3) 



TlLO q = 25,; 



They are displayed in fig. 4. In a finite spin cluster, how- 
ever, only the discrete waves with the g's restricted to the 
above values are possible. The spin-wave energies for a 
N = 6 system are indicated in fig. 4 by open symbols. For 
an antiferromagnetic 7V=6 ring, four almost degenerate 
spin waves are predicted, while for an antiferromagnetic 
N~6 chain the spin waves are organized into two groups, 
exactly as was just discussed. 

As pointed out already in the above, the rotational- 
band structure implies a classical spin structure in the L 
and E bands. This suggests to compare the results of 
various semi-classical theories (the details of our calcula- 
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Fig. 5: Calculated levels for an antiferromagnetic N — 6 spin- 
3/2 chain comparing the exact result and the results of the 
SWT calculations for J = -1.27 meV (LSWT = linear SWT, 
ISWT = interacting SWT, MLSWT = modified-linear SWT, 
MISWT = modified-interacting SWT, SLSWT = spin-level 
SWT) . The shaded areas highlight the two lowest L-band states 
and the _E-band triplets. 



tions will be reported elsewhere). The most widely rec- 
ognized semi-classical theory in magnetism is SWT. Here, 
one starts from the classical ground state (which in the 
present case is the Neel state with all spins oriented up or 
down) and accounts for quantum fluctuations by replacing 
the spin operators by Holstein-Primakoff or Dyson- Maleev 
bosons [13]. This yields a boson Hamiltonian TCb in real 
space, which is the appropriate starting point for finite 
spin systems as pointed out in Ref. [16]. 

In linear SWT, the boson Hamiltonian Hb is truncated 
to terms quadratic in the boson operators. The spin-wave 
energies are then obtained from a Bogoliubov diagonal- 
ization (in this work, the Bogoliubov diagonalization was 
performed numerically for all SWT calculations). Inter- 
acting SWT additionally includes the effects of 4th-order 
boson terms [13]. 

A serious drawback of these kind of theories is that 
they start from the assumption of an ordered, symmetry- 
broken ground state, which for finite clusters is obviously 
incorrect. As a result they are unable to reproduce the 
singlet-triplet gap, which is characteristic for finite spin 
systems. These theories hence miss the experimentally 
observed peak I. Modified SWTs have become popular 
to correct for this deficiency [16, 26]. Here, the rota- 
tional invariance is restored by imposing the conditions 
(Sf) = of zero on-site magnetization via Lagrange mul- 
tipliers. Modified-linear SWT goes up to quadratic boson 
terms, and modified-interacting SWT up to 4th-order bo- 
son terms. 

Finally, a further approach is suggested. For antifer- 
romagnetic rings (and similarly for other antiferromag- 
netic clusters such as Fe3o) the rotational-band struc- 
ture implies approximating the wavefunctions of the L 
and E bands by the spin levels IcxSaSbSAI), where Sa 
and Sb are the total spins on each of the two antifer- 
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romagnctic sublattices A and B (a denotes intermediate 
spin quantum numbers) [6, 8, 10]. The L band is ob- 
tained for Sa = Sb = SiN/2, while the space of the E 
band is spanned by the spin levels with Sa = SiN/2, 
Sb = SiN/2 — 1, and vice versa [10]. This suggests to 
diagonalize H in the space of the spin levels \qlSaSbSM) 
with the above values for Sa, Sb- This approach is called 
spin-level SWT. 

The results obtained for Ctq with all these techniques 
are shown in fig. 5 and compared to the exact values ob- 
tained numerically by exact diagonalization. Some gen- 
eral trends can be observed (which have been confirmed 
by studies on more spin clusters). First the E band is dis- 
cussed. For infinite spin systems the experience is that lin- 
ear SWT underestimates energies, while interacting SWT 
often yields very good results [13]. For Cr 6 a similar trend 
is observed. The linear SWT results are well below the 
exact energies, while interacting SWT comes rather close, 
within 8% accuracy for Cr 6 . The results of the modi- 
fied SWTs show a similar trend, i.e. modificd-interacting 
SWT yields a better agreement with the exact result than 
modificd-linear SWT, although not quite as good as inter- 
acting SWT. In contrast to these (bosonic) SWTs, spin- 
level SWT overestimates the spin- wave energies, by about 
30% for Cr 6 . As mentioned before, the modified SWTs 
yield a finite singlet-triplet gap Ao as does the spin-level 
SWT. A from modified-interacting SWT turns out 30% 
smaller than the exact result, while the spin- level SWT 
calculates it 50% too large (the performance of modified- 
linear SWT is even worse). 

Hence, we observe that the bosonic SWTs underesti- 
mate the energies while the spin-level SWT overestimates 
them. The interacting SWT yields the best results for the 
spin-wave energies, but is not able to reproduce a non- 
zero singlet-triplet gap (by construction). The modificd- 
interacting and spin-level SWTs both yield a non-zero 
value for the singlet-triplet gap A , but the overall ac- 
curacy is modest, with no clear advantage for one or 
the other. A conceptional advantage of spin-level SWT 
though is that it is rotationally invariant, and that it can 
be naturally extended to systems with, e.g., three sublat- 
tices such as Fe3o, where the non-collinear spin structures 
require special attention in the interacting SWTs. 

In conclusion, we have observed the standing spin-wave 
excitations in the molecular antifcrromagnctic finite-chain 
compound Cr 6 experimentally. We systematically investi- 
gated the capabilities of various SWTs in describing the 
elementary excitations of finite Heisenberg spin clusters. 
In our opinion, such investigations are of great value with 
regard to molecular nanomagnets where exact results are 
not accessible. The extendability of the rotational-band 
concept to spin systems such as finite antiferromagnetic 
chains underpins a certain generality. 
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